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Abstract 

We establish eigenfunctions estimates, in the semi-classical regime, for critical en- 
ergy levels associated to an isolated singularity. For Schrodinger operators, the 
asymptotic repartition of eigenvectors is the same as in the regular case, excepted 
in dimension 1 where a concentration at the critical point occurs. This principle ex- 
tends to pseudo-differential operators and the limit measure is the Liouville measure 
as long as the singularity remains integrable. 

Key words: Quantum chaos, Schrodinger operators, Equilibria in classical 
mechanics. 



1 Introduction. 



The problem we consider here concerns the asymptotic behavior of eigenvec- 
tors of a self-adjoint operator and follows the works of Colin de Verdiere [6] 
and Zelditch [13], on the basis of a result stated first by Schnirelman. We are 
more precisely interested in a proof of a microlocal concentration phenomena, 
in the semi-classical regime, as established in [1,7]. The adaptation to semi- 
classical analysis was done in [10], following a technic proposed by Voros [12]. 
We also mention [8] for a more general approach in the scattering setting. 

Consider a quantum operator P h , realized as a self-adjoint operator acting 
on a dense subset of L 2 (M n ). A typical example, studied in section 2, is the 
Schrodinger operator Ph = —h 2 A + V where the potential V is smooth and 
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bounded from below. If the spectrum of Ph is discrete in [E — e, E + e], a 
sufficient condition for this is given below, we can enumerate the eigenvalues 
in this interval as a sequence Xj(h) with finite multiplicities. We note ipj the 
corresponding normalized eigenvectors, i.e. 

P h rf = \j(h)^, ||i#||l" = i. 

Our objective is to establish eigenfunctions estimates: 



v J (a) = (Op w h (a)^i^ 
\j(h) -> E, 0+, 

where a G S' (lR 2ri ), so that by the Calderon-Vaillancourt Theorem: 
f^Op w h (a)f(x) = -^— J a £±l-,Z)ete-y*>f(y)dydZ, 

is bounded on L 2 (R"). Note that the statement of the problem is local w.r.t. 
E and we are here interested in the case of E — E c critical. Each Vj{a) 
measures the observable Op™ (a) in the state tpj. Let $ t = exp(tH p ) be the 
Hamiltonian flow of the principal symbol p of Ph- Interpreted as distributions, 
these measures are <3> t -invariant which easily follows from: 

e-^Op£(a)e*^,^) = (OpJ^e^ty* e^tyj) = uj(a). 



By Egorov's Theorem e h Ph Oph( a ) e h Ph is an operator of principal symbol 
(a o $ t ) and Vj is invariant under <3>j, up to 0(h). 

We recall that E is regular if Vp 7^ on the energy surface: 

and critical otherwise. When E c is critical, S Sc is not a smooth manifold. For 
E regular, T, E inherits a measure, invariant by $ t , given by: 

dz 

where dz is the Riemannian surface element. We note V(E) the associated 
volume of and we obtain a probability measure via: 

d» E (z) = ^ydLvol(z). 
Note that if E c is critical, dfx Ec has a sense if and only if 1 e L 1 (S Sc , dLvol). 
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Via a wave equation approach, substituting here the heat equation strategy 
of [6] on a compact manifold, the problem is related at the first order to the 
geometry of the energy surfaces. Accordingly, the integrability of dLvol has a 
strong effect on the asymptotic behavior of the sequence Vj. 

The micro-local concentration near a singularity was studied in [7] in dimen- 
sion one for a non-degenerate instable equilibrium where they proved that Uj 
converges to the Dirac mass at the equilibrium. We are first interested in the 
case of Schrodinger operators, but a generalization to pseudo-differential oper- 
ators provides more examples. Finally, we recall that in Riemannian geometry, 
e.g. for compact surfaces of negative curvature, the question to know if the 
full sequence converges to the invariant measure (quantum unique ergodicity) 
is still an open problem. 

Definitions. 

We define now the objects used in sections 2&3. If E c is a critical energy level, 
we pick an /i-dependant interval: 

1(h) = [E c -dh,E c + dh], d>0. (1) 

The associated counting function of eigenvalues is: 

T(ft)= E {^^) = W I \(h)el{h)}. (2) 

\j(h)€l(h) 

For A = Op^(a) a pseudodifferential operator of order zero, whose principal 
symbol is a e S^IR 2 ™) we put: 

T«(/0= E (A1>},il>})= E ";(«)■ (3) 
\j(h)ei(h) *j(h)ei(h) 

Observe that T(h) = Ti(/i) since for every quantization used in this contri- 
bution the symbol of the identity is I. 



2 Schrodinger operators. 



Let p(x,£) = £ 2 + V(x), where the potential V G C co (R n ) is real valued. To 
obtain a well defined spectral problem, we use: 

(Hi) There exists CeK such that: liminf V > C. 

oo 

By a classical result, Ph = —h 2 A + V(x) is essentially self-adjoint. Note that 
(Hi) is always satisfied if V goes to infinity at infinity. Let J = [Ei, E 2 ], with 
E 2 < liminf V. Since p _1 (J) is compact, the spectrum a(P h ) fl J is discrete 
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and consists in a sequence Xi(h) < \ 2 (h) < ... < Xj(h) of eigenvalues of finite 
multiplicities, if h is small enough. Next, we impose the singularity: 

(H 2 ) On T, Ec the symbol p has an isolated critical point z = (x ,0). This 
critical points can be degenerate but is associated to a local extremum of V 

V(x) = E c + V 2k (x) + 0(\\x - 4\\ 2k+1 ), keN*, 

where V 2 k, homogeneous of degree 2k, is definite positive or negative. 

The case k — 1, i.e. a non-degenerate singularity in dimension n, is treated in 
[1] without any extremum condition. In Fig.l the line in bold is the critical 
energy level attached to the top of a one dimensional symmetric degenerate 
double well. Observe the unstability near the recurrent critical point. 





Figure 1. Energies surfaces of V(x) = — x 4 + x 6 . 



To simplify notations we write z — (x,£) G T*IR n and zq for a critical point. 
The first result concerns the statistical behavior of the sequence z^(a). 



Theorem 1 Assume (Hi) and (H 2 ) satisfied. If n > 1 we have 

Hits J'.\ / : j (1(1 . 



r a (h) 

hZo+ T(h) 



but in dimension 1 we obtain: 

r r a (h) 



(5 ZQ ,a) = a(x ,0). 



These relations are statistical since when h — > the number of eigenvalues in 
1(h) tends to infinity, see Proposition 4. We define: 



K(h) = {j e N j Xj{h) e 1(h)}. 
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Following the results of [1,6,13,10], for n > 1 Theorem 1 implies that if <£>j is 
ergodic on S Ec there exists a density one subset L(h) C K(h) such that for 
all integer valued function h — > G we have: 

Km(^( fc) ,OpJf(a)^ (fc) >= / ad^. 
Here, density one simply means that: 

Since there is nothing new to prove, we refer to [1,10] for a precise study. More 
interesting, is the generalization of the result of [7]: 

Corollary 2 In dimension 1, assume that E^ c is connected. Then, for the 
weak * topology, we have Vj{h) ~ * d(zo), j(h) G K(h), uniformly as h — > 0. 

Hence, we obtain a concentration at z . Observe that 8{zq) is $ t -invariant. If 
a is simply a function, the quantum probability satisfies: 

lim Xj(h) = E C ^ lim \^ {h) \ 2 (x) = 5 XQ . 

The interpretation is as follows. In dimension one, the invariant measure on 
T,e c has a singularity: 

dz 

dLvol(z) ~ c-j-j rr, near z = (x , 0). 

M ^ — 1 1 

The measure is not integrable and the result has to be different. For n > 1, 
the singularity is integrable and an isolated critical point has no effect. This 
reinforce the universality of the Liouville measure in quantum ergodicity. But 
the case n — 1 is important since many problems, with symmetries, can be 
reduced to the study of such a singular Schrodinger equation. See e.g. [7] for 
an application in Riemannian-geometry. 

Preliminary remarks. 

The case of a local minimum of V is not really deep. Since £ 2 > 0, zo is a local 
extremum of p and is an isolated point of S^ c . According to the results of [3], 
the contribution of a minimum is significative only if n — 1. We consider now 
the non-trivial case of a local maximum of V, corresponding to an unstable 
equilibrium of the flow. Finally, since we use below the functional calculous, 
p has to be a symbol. But with (Hi) we can eventually modify the potential 
V outside of a compact, without modifying the main results. Hence no extra 
assumption is required. Similar comments apply for section 3. 

Proof of Theorem 1. 
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We use the semi-classical trace formula technic. This approach, analogous to 
the trace of the heat operator of [6], uses the propagator e~h Ph and a general- 
ization of the Poisson summation formula for this operator. Let G «S(M), a 
Schwartz function. To approximate T we define: 

>y(E e ,h,<p)= E <p( Xi{k) ~ Ee )- (4) 

\Xj(h)-E c \<e U 

This object can be treated by mean of Fourier integral operators (F.I.O.), see 
e.g. [9], and we refer to [1,2,3] for a detailed study of the trace. We recall that 
the Tauberian approximation concerns expressions: 

TS>) = E«i(^M^ z ^)- 

Under our assumptions, the behavior of T E h determines the behavior of the 
weighted counting functions: 

Nl d (h) = E «;(*)■ 

\Xj(h)-E\<dh 

Precisely, we will apply the results of section 6 of [1] for (Xj(h) = 1 or Vj(a). 
Strictly speaking, aj(h) > is required but, by a standard result of pseudod- 
ifferential calculus, we can modify the quantization to have Vj(a) > 0. This 
does not change the main results, see Eq.(14) below. 

To attain our objective, we can suppose that supp(</3) C [— M, M], M > 0. 
Let 6 G C^(]E C — e, E c + e[), such that = 1 in a neighborhood of E c and 
< < 1 on E. We localize the problem near E c by writing: 

l(E c , h, ip) = li(E c , h, if) + -f 2 (E c , h, if), 

with: 

,\j(h)-E c . 



l2 (E c , h,<p)= E &(Xj(h)M^- — j 

\\j(h)-E c \<e 

= Tr0(P,M^ T ^). 

Where the last equality holds by support considerations. By a classical result, 
see e.g. [2] Lemma 1, the term 7i = 7 — 72 satisfies : 

71 (£ c , h, if) = 0(h°°), zsh^O. (5) 

The Fourier inversion formula for 72 and the previous estimate provide: 

7 (£ c , h, ip) = ^Tr J e itJ t<p(t)eM-rfPhMPh)dt + 0(h°°). (6) 
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Next, with a function * e C^ (T*R™), with * = 1 near z , we write: 

72 (£ c , h, <p) = 7 Z0 (E C , h, <p) + 7rc g (^c, h, <p), 

where: 

lzo {E c ,h,ip) = i-Tr J e^p(t)^ w (x,hD x )exp(-^-tP h )Q(P h )dt, (7) 

and 7 rcg is simply the difference. The micro-local term 7 Z0 contains the con- 
tribution of the singularity and the discussion below determines if this term 
is dominant. For finitely many critical point on S^ c , we could repeat the pro- 
cedure. For the convenience of the reader, we recall the contributions of an 
equilibrium to the trace formula. 

Proposition 3 Assume (Hi), (H 2 ) and that p e C£°([-M,M]) ; M < M . 
If xq is a local maximum of V we have: 

lzo (E c ,h,p)~h- n +^ £ £ h^Mog(hrA j:l , m (<P)- 

m=0,l j,l£N 2 

If "^fc 1 ^ e N and n is odd then the top- order term is : 

C n , k \og(h)h- n+ % + % J \V 2k (ri)\-%dri J |t| n W-V(t)dt. 

§™-i M 

Otherwise the first non-zero coefficient are given by: 

h- n+ ^{T n ^p) J \V 2k (rj)\-^drj. 

gn-1 



This result is the contribution of an equilibrium since the distributional coef- 
ficients have a non-discrete support, contrary to the Weyl-term supported in 
t — and the contributions of closed orbits supported by the length spectrum. 
The distributions T n ^ and the universal constants C rit k ^ depend only on 
(n, k) and are explicitly determined in [4]. Mainly, we need the order w.r.t. h 
of these contributions determined by the functions: 

Wc (h) = h- n+ % + % log{h) j , j = or 1. 

First, we give a natural application of Proposition 3. 

Proposition 4 The microlocal counting function satisfies: 



T(h) 



2dV(E c )(2nh) 1 - n + o(h 1 - n ), ifn>l, 
HX[-d,d]) w c(h) + o(w c (h)), ifn = l. 
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Here A is the first non-zero distribution of Prop. 3 and X[-d,d\ the characteristic 
function of [—d, d, ] . 

Proof. By construction we have: 

lreg (E c ,h,<p) = i-Tr J e^0(t)(l - V w (x, hD x ))exp(-^-tP h )Q(P h )dt. 

K 

By the standard calculus on F.I.O. and an easy application of the stationary 
phase method, as h tends to we obtain: 

llcg (E c , h, <p) ~ (2 ^?Li Lvol (^c n su PP (i - *)) + o(h 2 - n )- 

Hence 7 reg always contributes at the order h 1 ~ n . 
Case of n > 1. We have w c (h) = o(h l ~ n ) and : 

Lvol(S Sc n supp(l - *)) < V(£ c ) < co, Vtf. 
It follows easily by shrinking the support of the cut-off \P that: 

7 (£ c , /i) = (2nh) 1 - n mV(E c ) + o(/i 1 -"). (8) 
Since the distributional factor is : 

0(0) = | </?(*)<&= (!,</?}, 
replacing ip by X[-d,d], y i a Theorem 6.3 of [1], provides: 
T(h) = 2dV{E c ){2nh) 1 ' n + o(h l ~ n ). 

Case of n — 1. Here the contribution of the critical point has a bigger order 
than the regular one. We obtain : 

-y(E c , </?, h) = w c (h)K(tp) + o(w c (/i)). (9) 

To apply the Tauberian argument of [1], we observe that our distribution 
A G iS'(R) can be represented by an element of L 1 1 oc (IR) and hence can be 
extended as a linear form on Cq(M) H L°°(M). We obtain: 



T(h) ~ A( X[ _ M ) Wc (h). (10) 

Which provides the desired result for n — 1. ■ 
As a matter of illustration, for n — k — 1, we have: 

logWl^o)|-/^ = ^^. (11) 

^ |^"(x )h 
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Which is the result established for T(h) in [1,7]. Observe that for n — 1, k — 1 
is the only case where a logarithm occurs and T(h) is slowly increasing w.r.t. 
k since for all k > 1 : 

r(h) ~ c/i^-5|y( 2fc )(x )r^. (12) 

Eigenfunctions estimates. 

We recall how to derive eigenfunctions estimates from the trace formula. First, 
to insert an observable A changes almost nothing. If II is the spectral projector 
on [E c — e, E + e], computing the trace in the basis ip^ and by cyclicity: 



' \^(h)-E c \<e a 



Since A is a bounded operator, if ip e <S(R) we can again smooth the problem 
via an energy cut-off B(P/ l ), with an error of order 0(h°°). Hence we can 
insert A = Op™ (a) in Eq.(7) and the results of Prop. 3 are the same after 
multiplication by a(z ). Similarly, the regular contribution changes via : 



0(0) 



a(z)(l - V(z))dLvol(z). 



By evaluation of the trace, we have : 

E K <p( Xj(h \~ E ° ) ~ c (ip)w(h)m(a) + o(w(h)), 

\\j(h)-E c \<e 11 

where w(h) changes only if n — 1. By Theorem 6.3 of [1] we obtain : 



Aj(fc)e/(/i) 
In particular this implies that: 

Substituting the correct expressions for these measures we obtain: 

- for n > 1 : m is a constant multiple of the Liouville measure. 

- for n = 1 : m is a multiple of the delta-Dirac distribution in z . 

Extraction of a subsequence. 

We chose a > and modify the quantization. Different choices are possible: 
Friedrichs quantization Op F as in [1] or the anti-Wick quantization Op AW as 
in [10]. These quantization are positive, i.e. 

a > (f,Op AW (a)f) > 0, V/ G C °°(R"). 
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Since Op^(a) — Op h (a) is of order -1, we obtain : 

(Opjfta) - Op£ w (a))tf) = 0(h), (14) 

and we can work with this positive operator. For n > 1, under the condition 
that $ f is ergodic on S^ c , the extraction of a convergent subsequence of density 
one is the same as in [1,6,10] to which we refer for a detailed proof. For n — 1, 
if E £ c is connected, there is only one probability measure invariant by $ 4 and 
the full sequence converges to 5 ZQ . Once the result is established for a positive 
symbol it can be extended by linearity to any a e S' (1R 2 ™). 



3 Pseudo-differential operators. 

The case of pseudo-differential operators provides more explicit examples. Let 
Ph = Op™(p(x,£,), obtained by Weyl quantization, where the symbol p is real 
valued and smooth on T*IR n . In general position, one can also consider h- 
dependent symbols J2 h j Pj(x,^), see [10]. But, to simplify, we consider only 
the homogeneous case. As above we impose: 

(Al) There exists e > such that p _1 ([i? c — e , E c + e ]) is compact. 

As in section 2, a(P h ) n [E c — e, E c + e] is discrete. A fortiori, (A\) insures 
that T,e c is compact. Next, we chose an homogeneous singularity : 

(A2) On T,e c , p has a unique critical point z = (xo,£o) an< ^ near z o '■ 

p(z) =E C + p k (z) + 0(\\(z - z )\\ k+1 ), k > 2, 

where p& is homogeneous of degree k w.r.t. z — z . 

Strictly speaking, one could consider k = 2. But this case is precisely treated 
in [1]. The case of a critical point which is not an extremum is technical 
because the singularity is transferred on the blow up of zq. To obtain a problem 
that can be explicitly solved, we consider the following hypothesis inspired by 
Hormander's real principal condition: 

(^3) We have Vp k ^ on the set C{p k ) = {9 e S 2 "" 1 / p fc (0) = 0}. 

For example ps(x,£) = x 3 — £ 3 is admissible and p(x,£) = x 3 — £ 3 + x 4 + £ 4 
satisfies all our hypothesis for E c = 0. 

Remark 5 With (.A3), contrary to the case of a local extremum, zq is not 
an isolated point of which imposes to study the classical dynamic in a 
neighborhood of z Q . The study of singularities like in (A3) is detailed in [11] 
chapter 4 to which we refer concerning the integrability of dLvol. 
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As in section 2, it is sufficient to study the local problem 7 2o defined in Eq.(7). 
The contributions to the trace formula are : 

Proposition 6 Under (Ai) to (A3), we have an asymptotic expansion: 

00 

lzo (E c ,y>,h)~h%- n £ Y.h^og(h) m K hm ( V ), 

m=0,l j=0 

where the logarithms only occur when (2n + j)/lc 6 N* and Aj >m e «S'(R). 
As concerns the leading term we obtain: 

(1) If k > 2n (non-integrable singularity on T,e c ) we have : 

j Z0 (E c ,<p,h) ~ h^- n Ao, (ip) + 0(h^- n log(h)), ash^O, 
where Ao,o is a universal distribution. 

(2) If the ratio 2n/k G N we obtain logarithmic contributions : 

lzo (E c ,<p,h) ~ /^- n log(/i)A ,i(y?) + 0(hT- n ), ash^O, 

(3) For 2n > k, 2n/k ^ N the result is as in (1) with a different distribution. 



These results describe very precisely the singularity at z . But this is not our 
purpose here and we refer to [5] for a detailed formulation of these contri- 
butions. For n — 1, k — 2, the case (2) agrees with section 2 and allows to 
recover some results established in [1,7]. 

Application to microlocal measures. 

The proof is exactly the same as in section 2. The main difference is that 
the singularity on Se c can be of arbitrary order. In our setting, according to 
Prop. 6 the top order coefficient changes if and only if we have: 

2n 1 2n 1 

n < 1 - n — < 1. (15) 

k k 

If k < 2n the singularity is integrable and contributes at a lower order com- 
pared to h 1 ^ n V(E c ). But if k > 2n, which corresponds to a non-integrable 
singularity for d/i Ec , the main term changes. To summarize, we obtain: 

/ adfi Ec , for k < 2n, 
a(zo), for k > 2n. 

Contrary to section 2, observe that for k > 2n and if n > 1 we do not obtain 
the convergence of the full sequence Vj(h), j(h) € K(h), to the dirac-mass at 
the equilibrium. The obstruction is that an invariant probability measure can 
be supported by the closed orbits of S Sc . 

Comments. 



= < 
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From these 2 families of example the conclusion is that the limiting measure 
changes only if £e c carries a measure such that 1 ^ L\ oc {Yie c , dLvol). Inter- 
preted as a quantum measurement, one obtain a very precise localization : if 
a = around z the limit is the Liouville-measure but if a(z ) ^ the limit 
strongly differs. 

An interesting problem would be to study the repartition in presence of 2 equi- 
libria zi, z 2 on T lEc of the same nature and with a non-integrable singularity. 
In this case any convex combination : 



provides an invariant probability measure. A natural question is to determine 
if the limiting measures are equally distributed between z\ and z 2 . 
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